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Generalized Elastic Model: thermal vs non-thermal initial condi 
tions. Universal scaling, roughening, ageing and ergodicity. 



(N 

o 



^ 









o 

(N 
> 

0^ 
(N 
en 

rn 

O 



X 



Alessandro Taloni^'^ , Aleksei Chechkin^'^ and Joseph Klafter^ 

^ School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel 

2 CNR-IENI, Via R. Cozzi 53, 20125 Milano, Italy 

^ Akhiezer Institute for Theoretical Physics, NSC KIPT, Kharkov 61108, Ukraine 



PACS 05 . 40 . -a - Fluctuation phenomena, random processes, noise, and Brownian motion 
PACS 02 . 50 . Ey - Stochastic processes 
PACS 87 . 10 . Mn - Stochastic modehng 

Abstract -We study correlation properties of the generalized elastic model which accounts for 
the dynamics of polymers, membranes, surfaces and fluctuating interfaces, among others. We 
develop a theoretical framework which leads to the emergence of universal scaling laws for sys- 
tems starting from thermal (equilibrium) or non-thermal (non-equilibrium) initial conditions. Our 
analysis incorporates and broadens previous results such as observables' double scaling regimes, 
(super)roughening and anomalous diffusion, and furnishes a new scaling behavior for correlation 
functions at small times (long distances). We discuss ageing and ergodic properties of the gen- 
eralized elastic model in non-equilibrium conditions, providing a comparison with the situation 
occurring in continuous time random walk. Our analysis also allows to assess which observable is 
able to distinguish whether the system is in or far from equilibrium conditions in an experimental 
set-up. 



Introduction. — Continuum linear systems enjoy 
an evergreen popularity among the scientific community. 
This is due to the simplicity of their formulation, to their 
solvability and to their apparent capability of capturing 
and reproducing the dynamics of complex natural phe- 
nomena. This is the case, for instance, of polymers [1], 
membranes [21l3], interfaces [IHl], surfaces [8|, single-file 
models [9] , all systems that are governed by the following 
stochastic differential equation 



d_. 

dt' 



h{x,t) = / dVA(x-f') 



a^ 



d\ 



-h{x',t) + T]{x,t), 



(1) 

which has been introduced as to generalized elastic model 
(GEM) [TU] . The motion of the D-dimensional stochas- 
tic field h {x, t) on the d-dimensional substrate x is 
driven by the hydrodynamic interactions, embodied by 



A(r) 



z/2 



by the fractional Laplacian 



0" 

d\s\- 



— (— V^)"'" [T31, and by the thermal noise random source 
r){x,t) fulfilling the fluctuation-dissipation relation (FD), 
i.e. {rjji£,t)r]kix',t'))^2kBTA{x-x')SjkSit-t'). We 
hereby consider only systems satisfying z > d, for which 
the interface is called rough [6]. Moreover, systems with 



local (or screened) hydrodynamic interactions are charac- 
terized hy A(r) = S (r) . 

However, the system's dynamics and the ensuing macro- 
scopic observables strongly depends on the initial condi- 
tions of the equation ([T]). For instance, assuming the sys- 
tem in a stationary state at t = corresponds to take a 
polymer in a coiled relaxed configuration before starting 
any measurement. The same can be thought for a floating 
membrane or a surface, or a rough interface, which has 
reached the thermal equilibrium condition much before 
than any experimental observation led off. We will refer 
to these situations as thermal initial conditions. On the 
other hand, the scientific interest is often directed towards 
the relaxational properties of the system under study: a 
common self-question which an experimentalist (or a the- 
orist) asks is "how do i know/characterize the behavior 
of a system which is relaxing to its equilibrium configura- 
tion?" . Think, for instance, to the polymer relaxational 
dynamics after the translocation through a nanochannel 
or a nanoslit, to the evolution of a membrane from a flat 
initial condition, to the growth of an interface on a flat 
substrate, to the diffusional properties of an assembly of 
1-dimensional Brownian particles equally spaced at t = 0. 
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These initial conditions go under the name of non-thermal 
initial conditions and are characterized by h (if, 0) =0 Va; 
[To]. Moreover non-thermal fluctuations, often referred 
to as "cytoskclctal fluctuations", have been shown to be 
responsible for the non-equilibrium dynamics in many bi- 
ological systems, such as traveling waves due to protein 
activity in a flexible membrane, as well as in micropipet ex- 
periments using an activated membrane surface (see [11] 
and references therein) . The aim of this letter is to show 
how the scaling properties of a generic physical observable 
depends on the chosen initial condition, whether thermal 
or non-thermal. 

Mathematically the question is well-posed. In the case 
of thermal initial conditions we assume that the system 
reached the stationarity at t = — oo: the natural conse- 
quence is the use of Fourier transform in time and space, 
i.e. h(g,w) = /+^d'^a;/+^dih(f,i) e-*(«-^-"*). 
The Fourier-Fourier transform of the ther- 
mal noise indeed reads {rjj {q,uj)i]k{(f ,io')) = 

cf)S{uj + Lu'), where 
< a < d. In- 
stead, in case of non-thermal initial conditions we 
will make use of the Laplace transform in time: 
h{q,s) = /+^d''x/o^°°dth(x,t) e-*(«-^)-^*, then 
the noise Fourier-Laplace transform is given by 

The local hydrodynamic situation is achieved by setting 
formally A = const and a = d in the corresponding 
long-ranged hydrodynamic expressions; however, this 
substitution is not to be intended as a limit [TUlfT^ . 
The results of our analysis lead to the emergence of an 
universal scaling framework, and ensuing time scale, char- 
acterizing the observables' behaviours in both thermal 
and non-thermal initial conditions. Within this context, 
we analyze different scaling regimes displayed by the 
correlation functions for systems presenting long range or 
local hydrodynamic interactions. Furthermore, we discuss 
the ageing and ergodic properties of a system during its 
non-equilibrium relaxational phase. 

Tv^ro-point two-time h-correlation function. 

We hereby derive the following two-point two-time corre- 
lation function for the stochastic flcld component hj (x, t) 



{[h {x, t)~h (x, 0)] [h (f , t') ~ h (f, 0)]) 
{5th{x,t)5t.h{x',t')), 



(2) 



where we implicitly dropped the index j . We aim at fur- 
nishing its scaling expression for systems characterized by 
long range or local hydrodynamic interactions, whether or 
not the initial conditions arc taken to be stationary. 

Thermal initial conditions. Our starting point is the 
solution of the GEM ([T]) in its Fourier representation 
for thermal initial conditions, which reads h{q,uj) = 

■q {q, Lu) I { —iu + A |q|"^"~ ) . Hence, using the Fourier 
transform of the noise correlation function we derive 



{h (q, u) h {<f,Lu'))th = 2kBT(27r)^+^^is\ 
S{q'+q')6{uj + u}'), 



uj'^+A'^\q\'> 



(3) 



where 7 — 2{z -\- a — d) . This expression will play a central 
role in the analysis that will follow. After applying the 
inverse Fourier transform in space and time, we get 



•^-^ If-ff-'^/'x 



{h {x, t) h (f , t'))th - (.3^)^/2 

/- d \q\ \qf''-'- J,/,_i M \X-X'\) e-^\^'''V-''\ 

(4) 
where Jd/2-i{f) represents the Besscl function of order 
d/2 — 1. Inserting ^ into the definition ([2]), the two- 
point two-time correlation function can be arranged in the 
following general compact form 



{5th{x,t)5t,h{x' ,t'))th ^ ji^ 



f 



""^^ ^x-xT^'x 



t 


+ f 


[r{\S-S'\)\ 


-f 


t{\S-S'\) 


[t{\S~S'\)\ 



(5) 



where 



T[ X - X ^ 



is defined as the correlation 



time of the distance |x — x'j [12]. The scaling function 
/ [u] is given by 

/•oo 

f[u]^J dyy'^'-^-Jd/2-i{y){l-e-y'"'^) (6) 

and will be the subject of the upcoming investigation. As 
a matter of fact, the analysis of (jS]) reveals two different 
scaling behaviors whether u ^ 1 or u^ 1. Analyzing the 
scaling of the correlation function is not a mere mathe- 
matical exercise, since in real experiments the scattering 
data of polymers, membranes or rough surfaces require an 
expression of cq.Q at small and large times or respec- 
tively, for long and short distances J51 [T^[TB1[T7] . To show 
this, we first manipulate (jG]) to get 



/N 



du'u 



dvv''-'"^Jd/2-i(- 



-'2/7 



(7) 

For u ^ 1, systems characterized by long range hydrody- 
namic interactions behave unlike systems where hydrody- 
namic interactions are purely local. In the former case, we 
can safely replace the exponential in ([7]) by 1, achieving 



■'LI -p / d—a 



(8) 



On the other hand, when the hydrodynamic interactions 
are completely screened, the exponential in ([7]) must be 
expanded to the first order. The ensuing integral can be 
computed yielding as a result 



/H.2— i.„(f)r(£)r 



z + d\ 2 /o\ 
u , (9) 
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Title 



when z ^ 2m {m G N). If z = 2to, the scahng function is 
exponentiahy smaU so that (up to constant factors) 



/ [u] ex u^'+ig-" 



(10) 



this resuh extends to a general system the cases 2; = 2,4 
and d ~ 1 studied in [14]. The novelty of eq.® and 
([9]) and their experimental relevance will be highlighted 
in the next section. In the opposite limit u 3> 1, there 
is no difference among systems with long range or local 
hydrodynamic interactions: this can be seen by expanding 
to the first order the Bessel function and calculating the 
ensuing integral, the result reads 



/M 



2i-rf/2r(i-/3) 



r(i) 



(11) 



where ^ = 2(z - d)/7, < /3 < 1 [TU|. We notice that, 
thanks to the definition of the correlation time t {\x — x'\), 
the long time limit expression (|lip entails a cancella- 
tion of the spatial dependence of the two-point two-time 
correlation function ([5]), i.e. {5th{x,t) 5t'h{x' ,t'))th — 
{5th (x, t) Sfh {x, t'))th, where 



{5th{x,t)5t,h(x,t'))th^K tP+t'P-\t-t'f (12) 



and K 



2kBTTr'''^ A^r{l-/3) 



as already pointed out in 



(27r)''r{d/2) z-d 

[6IIIO). Therefore the characteristic time r(|a; — x'|) can 
be seen as the time up to which the dynamics of two 
distinct probes in x and x' is uncorrelated, indeed for 
t ^ t{\x — x'\) the autocorrelation function (fT2|) coin- 
cides with the two-point two-time correlation function 
([5]). Alternatively, one can define the correlation length 
^(i) = (Ai)2/^ [5H7l[ni[16H22] and say that two probes 
are correlated when ^(i) exceeds the distance \x — x'\. 

Non-thermal initial conditions. When the GEM ([T]) 
starts from non-thermal initial conditions, its solution 
in the Fourier-Laplace space is given by h ((f, s) = 



T](M,s)l(s^A\q\ 



z-\-a — d 



from which it follows 



{hiq,s)hi<f,s'))nth^2kBT{2Tr) 



{s+A\q\-'){s'+A\q\-'}^ 



d ^(?+r) 

s-\-s' 



(13) 



thanks to expression for the Laplace transform of the noise 
correlation function. Making an inverse Fourier-Laplace 
transform in space and time yields 



{Sth{x,t)5t'h{x',t'))nth 



'^^^ ^x-xT^x 



I 



_t±t' 



t{\x — x'\) 



f 



(2^) 



T ( I X — X' I ) 



(14) 



Thus we trace the previous analysis leading to the ex- 
pressions dll), ®, (Uni) and dll]), in the limit u < 1 
(<±<' < t{\x-x'\)) and m > 1 (i±t' > T(|f-.f'|)) 
respectively. In particular {5th{x,t)5t'h{x\t'))nth ^ 



{Sth{x,t)St'h{x,t'))nth for i ± i' >• t {\x — x'\), where 



{5th [x, t) St'h (x, t'))ruh = K[{t + t'f -\t- tf ] (15) 

As an example, let us discuss the situation of local 
hydrodynamic interactions (A(|7^) =(5(|7^)) and z = 2, 
d = 1 and A = 1 in ([1]). This corresponds to the equa- 
tion for Rouse polymers [T], for the Edward- Wilkinson 
chains [4] , for the attachment-detachment diffusion model 
of fluctuating interfaces [5], for single-file systems [5], 
and it is also known as ID diffusion- noise equation |15] . 
In this case we can compute the expressions ^ exactly: 

f[u\ = \/2ue~'^ — V^ ^''/^ ( V 51 ) ' where erfc denotes 
the complementary error function. Expanding erfc for 
small and large arguments gives f[u] ~ {2u)'^^^e~^ if 
M ^ 1, and f[u] ~ v2m if u ^ 1, which corresponds to 
eq. (fTO| and eq.l^ respectively [9l[T4l[T6|. 

Two-point one-time correlation function. We 

now turn to the analysis of the following correlation func- 
tion 



{[h (x, t)-h (0, t)] [h (f, t)-h (0, t)] 
{S^h{x,t)S^'h{x',t)). 



(16) 



Our intent is to show the space scaling properties of the 
GEM (P at a given time t. 

Thermal initial conditions. If the entire system is in 
a stationary state, from (j3]) we can derive the structure 
factor or power spectrum J18j : 



{h{q,t)h{(f,t))th^kBT{27T 



\qr ' 



(17) 



Making an inverse Fourier trasform in space, the previous 
expression gives 



{h (x, t) h (x', t))th = ksT 



(2T)d/2 

Jd/2-1 {\q\ l^-^'l), 



(18) 



which is also obtainable from eq.Q by setting t ~t' . This 
allows to recast the expression ([T6l) as follows 



{6xh{x,t)Sx,h{x',t))th 



(19) 



(20) 



ji^ {tth [x\ + ith [s'] - ith [x - .f']} . 

The analytical form of the function £th [x] is 

^th [x] = f^ \x\'-'^J^dyy'^-'~^x 
{l-r(d/2)(|)'-''/'j,/2_i(2/)}. 

The analysis of ([20]) reveals that the integral in y is con- 
vergent if the condition z < d -I- 2 is satisfied. As a con- 
sequence, systems satisfying the condition z < d + 2, are 
characterized by a correlation function p^ which is that 
of a fractional Brownian motion, with the time replaced by 
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the spatial coordmate \x\, and the Hurst exponent H^ = 
{z — d)/2, as firstly reported in [7]. On the other hand, if 
z > d + 2, the divergence for y — >• requires the introduc- 
tion of a lower cut-off ^ , where L represents the maxi- 
mum size of the system. This entails £th [x] oc L^^''^^ \x\ 
which, plugged into ea.([T9l). points out that the route 

mean square difference y {[h{x,t) — h{x',t)] )th is pro- where 
portional to the distance |a?— a;'| for any system for which 
z > d+2, as already known for bending- energy dominated 
membranes [3] and for fluctuating interfaces [7]. 
Now, recalling the definition of the correlation time r (|x|), 
we can summarize the obtained results by recasting the 
two-point one-time correlation function (J19p in the same 
fashion as the two-point two-time correlation function ([5|) 

This can be done by expressing ^t/i [x] = '-^'^^ "■■ '^ 
with 5t/i[u] = gth ~ const ii z < d 



A gth ^YiWV) 
2, and gth[u] oc 
y2(z-d-2)/7 jf ^ -^ d + 2. Thus, the morphological tran- 
sition occurring atz = d-|-2 [51[7], turns out to be a 
general property satisfied also by systems with long range 
hydrodynamics such as membranes, polymers or viscoclas- 
tic surfaces. In analogy with fluctuating interfaces we will 
refer to the systems for which z < rf -I- 2 as Family- Vicsek 
systems, and to the systems fulfilling z > d + 2 as super- 
rough systems [T51ITO] . 

Non-thermal initial conditions. In case of non- 
thermal initial conditions, the correlation function 
{h {x, t) h [x' , t))nth can be obtained from the power spec- 
trum 



In analogy with the thermal case, the two-point one-time 
correlation function (IT51) can be casted as 



feflT 



{6^,h {x,t) Sx>h {x' ,t))nth = (2^)d/2 
inth [x', t] - (nth [x-x',t]} 



{I 



7ith 



x, t] + 



(23) 



d-2-l. 



(nth [x, t] = ^fjj^ Ix]"" J^ dy y' 

{i-r(d/2)(|)'-''/V,/2_i(y)}( 



1 



„-l/-^2t 



(24) 

2/7 



After the following change of variable v — y ( .^ngi\ ) 
the short and long time analysis of eqs.([Ml) and 1^^ can 
be performed. For short times, t <^ t (|a?|), to the leading 
order we get 



{6,h{x,t)S,,h{x\t)),uh-K{2t)^, 



(25) 



implying the cancellation of the spatial dependence in 
((23| . Eg. ([25)1 is exactly the mean square displacement 
of a tracer (probe particle) in a generic point x, when sys- 
tem starts from non-thermal initial conditions |10| . A 
different situation occurs whenever one is concerned with 



the long time limit, t 3> t (|x 



(nth [X,t] = |.f| 

g„th M as 



'.-d 



gnth 



t(|5i: 



), of ([M|) . First, we express 
. then we cast the function 



{h{q,t)h{q',t))„,, = kBT{2nr^-^^±fi (l - g-^^'^l** 

(21) 
which follows from eg. (1131) . Making an inverse Fourier 
transform in space we get 



gnth M 



<Jth 



- const X u 

2 z~d~2 

U ~i 



1- 



z-d-2 



z <d + 2 

z>d + 2 



(26) 



{h{x,t)h{x',t))nth 



knT 



f 



2t 



(27r)'*/2 ' ' lT[\x 

(22) 
due to the definition (j6]). Alternatively, the expression 
((221) follows immediately from cg. p^ setting t = t' . Any- 
how, the quantity (j22p is of experimental importance for 
the interpretation of scattering data, indeed it is straight- 
forwardly connected to the dynamic structure factor. By 
instance, in case of rough surfaces [HUH], a widely used 
form of (|22p for short times (long distances) is exponen- 
tial, in agreement with the expression (|10p [T7]. How- 
ever, we point out that this behaviour represents only a 
marginal situation (z = 2m) for systems characterized by 
local hydrodynamic interactions. As a matter of fact, the 
expression ([9]) yields a more general result: the long dis- 
tance behavior of ([2^ . for any z ^ 2m, exhibites a de- 
cay oc \x — x'\~^~ . On the other hand, eq.(j8|) furnishes 
a testable new prediction for the scattering data coming 
from long ranged hydrodynamics viscoclastic systems. 



This quantity is intimately connected to the emergence of 
the so-called anomalous scaling or anomalous roughening 
[51 [T51ll9j . This arises whenever the local width w'^ {l,t) ~ 
{{[h{x,t) — {h{x,t))id] );d) scales differently from the 
\ global width W'^ {L,t) = {{[h[x,t) - {h[x,t))Ldf)Ld). 
For t > t{L) the global width reaches the value 
W^ {L,t) oc L?^, where x = ^^ is the roughness expo- 
nent [22] , which coincides with the spatial Hurst exponent 
Hx [21]. On the other hand, the local width scaling law 
is given by w^ [l,t) - l^^'-- for t (0 < i < r (L). The 
local roughness exponent xioc can be, in general, different 
from X'- whenever it happens, a system is said to present 
anomalous roughening. The anomalous scaling has been 
detected in a large collection of experiments and models 
(see [I8H2T] and reference therein). In our framework 
it stems from the scaling properties of the function gnth. 



indeed 



{l,t) 



F 



gnth 



r(l) 



which gives xio 



X 



for Family- Vicsek systems and xioc = 1 for super-rough 
systems [I^. Our analysis extends this result to long 
range hydrodynamics systems, showing that anomalous 
roughening is a general phenomenon not only restricted to 
the domain of rough surfaces. By instance, flexible Zimm 
polymers [l] (z = 2, a = 1/2, d = 1 in cq.([T])) arc believed 
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to fall in the Family- Vicsek universality class, while fluid 
membranes and semiflexible polymers [11|3] (z = 4, a = 1 
and d = 2 and 1, respectively) should exhibit anomalous 
superrough scaling. As a matter of fact, the correlation 
function (|16p has been studied in [5] for membranes start- 
ing from thermal initial conditions, but its behavior has 
never been addressed in case of non-thermal initial condi- 
tions. 

Time vs ensemble average. — Our analysis now 
turns to the effects that the specific initial conditions may 
have on the time average of an observable, which is func- 
tion of a probe particle trajectory only, i.e. of a tracer 
particle placed at position x. In general, it is of a very 
broad interest to assess the correctness of the ergodic hy- 
pothesis in real systems, i.e. if the time and ensemble 
average coincide. In the specific, it might be extremely 
useful for an experimentalist to establish whether the time 
average of a one-point observable can be representative of 
the thermodynamic state of the entire system, if in or out 
of equilibrium. Firstly, we define the observable to be the 
square displacement Sfh{x,t) = [h{x,t) — h{x,Q)] . We 
then define the time average over a trajectory of length T, 
as 



■r 



5|/i(x,i) = 



T-t 



T-t 



dto [h (x, t + to) - h {x, to)Y 



(27) 

In general we expect that (j27|) is a fiuctuating quantity 

dependent on the chosen stochastic trajectory: it is then 

j- 

legitimate to take its ensemble average {6fh(x,t) ) and 

compare it with {S^h{x,t)) [5^. 

Thermal initial conditions. If the system is in ther- 
mal equilibrium, the probe dynamics is ergodic, indeed 
the stochastic motion of a probe particle in x is governed 
by a fractional Langevin equation with stationary frac- 
tional Gaussian noise [10], then limj-^oo S'^h{x,t) 
{Sfh{x,t))th = 2KtP {FigWia)). 



-T 



Non-thermal initial conditions. 



In order to calculate 



the ensemble average of p7|) . i.e. {6th {x,t) )nth, we 
need to get the probe MSD between times t + to and to: 
{[h{x,t + to) — h{x,to)] )nth- From (|15p it is immediate 
to obtain 



{[h {x, t + to) - h {x, to)] )nth = 

K{[2{t + to)f + {2toY + 2tP - 2{t + 2to)^} 



(28) 



At first, we notice that the MSD exhibits ageing since it 
displays a strong dependence on t and tg when both are 
large [15]. We then study its behavior for to '^ t and 
to >i 



{[h (f , t + to)-h {x, to)f)nth 



K{2tY 
2KtP 



to-^t, 
to > t. 



Therefore, for i <C to <C t{L), where t{L) represents the 
characteristic relaxational time of a system of size L, the 
probe's diffusion is like that for thermal initial conditions, 
albeit the entire system is in a non-stationary regime. In 
this condition an experimentalist could not decide whether 
the system is in equilibrium or not. These results are 
supported by numerical simulations of a single file sys- 
tem (see Fig[TJb)). An illustrative description of the ob- 
tained results may come from the comparison of ([28]) with 
the corresponding quantity of a subdiffusivc continuous 
time random walk (CTRW) process ^{t). It is known that 
CTRW MSD exhibits ageing P^H^. In the same limits of 



eq. (f29l) one has {[({t + to) - C{to)]) ^ f" for to < t and 

{[^{t + to) - ^(to)f) - to^^t for to > t, where ^i is the 
exponent of the waiting time distribution ?/'(t) ^ t~^~^ 
{0 < fi < 1). Thus, the MSD grows subdiffusively in the 
limit to/t — > but normally when to/t — > oo, which is 
at odd with the corresponding limits in ([29]) . where the 
difference is only in the prefactor. Moreover the ageing 
properties of the MSD in CTRW are very different from 
those of (P5|) : in particular, in the regime to ^ t, diffusion 
is slowed down in CTRW as tg increases [25], while the 
probe diffusion in case of non-thermal initial conditions 
loses the dependence on to, to the main order. 
We now turn to the ensemble average of (pT)). Integrating 
eg.ipS]) one gets 



tT, 



{dfh{x,t)' )nth = Kx 

2/^(r'^+i+(r-t)'^+i-t'^+i)-(2r-t)''+i+t''+i 
(i+/3)(r-t) 



2^3 



(30) 



from which it is apparent the dependence on the trajec- 
tory length T. This is indeed shown by the outcome of 
the numerical simulations of single file system, perfectly 
reproduced by the theoretical prediction (PU]) (Fi^ljc)). 
The natural requirement in experiments is to take the 
limit t <^ T- We then expand expression (PU)) for small 

values of the parameter t/T obtaining {5fh {x,t) )nth — 



2Kt(^ ( 1 



-1 t 



Therefore, if the system is pre- 



2(1-1-/3) T J 

pared in non-thermal initial conditions, in the limit of 
t/T -^ 0, the ensemble average of ((27)) tends to the value of 

the ensemble (or time) averaged MSD at equilibrium, i.e. 

f 

(Sfh {x, t) )nth -^ (Sfh {x, t))th- Thus, it would be impos- 
sible for an experimentalist to assess whether the system is 
in equilibrum or far from it, just by looking at the probe 
time averaged MSD when the condition t <C T is ful- 
filled. However, when t approaches T, the T-dependence 
appearing in pop would become apparent and it would 
constitute the signature of the non-equilibrium thermody- 
namical state of the entire system (Fi^ljc)). Finally, the 

expression ([50]) is in contrast to that obtained in CTRW 

■r. 



(29) 



model, for which {5fS,{t) ) - T'^^H when t-^T ^. 

Conclusions. — In this letter we established a univer- 
sal theoretical framework for the correlation functions in 
a generalized elastic model. Within this framework, the 
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Fig. 1: (Color online) MSD of a tagged tracer in a single-file 
system {z — 2, d — 1 and A = 1/^, with ( = 1.0 the damp- 
ing [9]). Simulations are performed integrating the stochas- 
tic dynamics of a file of A'^ = 2000 Brownian {ksT = 1.0) 
hard-core pointlike particles, moving along a ring of length 
L = 20000. Thermal initial conditions correspond to pick a 
random file's configuration from an uniform distribution at 
f = 0, whereas non-thermal initial conditions correspond to 
take particle equally spaced of a distance L/N {L/N = 10 in 
the simulation). Panel (a): MSD ensemble average for ther- 
mal (solid upper black line) and non-thermal initial conditions 
(solid bottom red line), plotted against the theoretical predic- 
tions 2Kt^^ and K{2t)^ respectively (dashed lines). Panel (b): 
{[h {x, t + to) — h {x, to)] )nth for different to (solid lines), upper 
black dashed line corresponds to 2Kt^ , bottom red dashed line 
to K{2t)^ , while middle color dashed lines stand for ea. H28[) . 

7" 

Panel (c): {Sfh{x,t) )nth for different trajectories's length 
T. Solid lines show the outcome of the numerical simula- 
tions, dashed color lines represent the corresponding theoreti- 
cal ea. (|30[) . Upper and bottom dashed lines are the same as in 
panel (b). Simulation results are averaged over 40000 different 
probe's trajectories. 



scale invariance and the scaling regimes attained by any 
physical observable emerge naturally, getting a physical 
meaning in terms of the correlation time t. On one hand, 
we found new scaling behaviors of the correlation functions 
at small times (large distances) providing testable predic- 
tions for the behavior of the dynamic structure factor in 
case of rough surfaces and viscoelastic systems. On the 
other, we showed that anomalous roughening is a phys- 
ical phenomenon that can be detected also in processes 
characterized by long-ranged hydrodynamic interactions. 
Moreover, we analyzed ageing and ergodic properties of 
the generalized elastic model performing time and ensem- 
ble average of the squared displacement. Our results point 
out that the probe's MSD attains its equilibrium value al- 
though the entire systems is far from it. 
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